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PREFACE 


This  report  summarizes  the  work  completed  during  the  third  year  of  a  project  focused  on  im¬ 
proving  the  WBMOD  model  of  ionospheric  scintillation.  This  work  is  part  of  a  larger  effort 
which  has  the  overall  objective  of  providing  the  USAF  Air  Weather  Service  with  an  im¬ 
proved  capability  for  providing  support  to  their  customers  who  are  impacted  by  ionospheric 
scintillation  effects. 

We  would  like  to  thank  Dr.  Santimay  Basu  of  the  USAF  Phillips  Laboratory  (PL/GPIA) 
for  his  guidance  and  suggestions  for  improving  the  model,  Dr.  Fred  Rich  (PL/GPSG)  for 
providing  the  DMSP  SSJ/4  particle  boundary  data,  and  Dr.  Ed  Fremouw  (NWRA)  for  many 
useful  suggestions  and  numerous  discussions  over  the  course  of  this  work. 
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1.  Introduction 


Many  modern  military  systems  used  for  communications,  command  and  control,  naviga¬ 
tion,  and  surveillance  depend  on  reliable  and  relatively  noise-free  transmission  of  radiowave 
signals  through  the  earth’s  ionosphere.  Small-scale  irregularities  in  the  ionospheric  density 
can  cause  severe  distortion,  known  as  radiowave  scintillation,  of  both  the  amplitude  and 
phase  of  these  signals.  A  basic  tool  used  in  estimating  these  effects  on  systems  is  a  computer 
program,  WBMOD,  based  on  a  single-scatter  phase-screen  propagation  model  and  a  number 
of  empirical  models  of  the  global  morphology  of  ionospheric  density  irregularities.  A  recent 
validation  of  the  WBMOD  model  showed  that  it  was  deficient  in  a  number  of  areas.  Much 
of  this  deficiency  could  be  traced  to  the  limited  coverage  of  the  data  used  in  developing 
the  model  (a  single  station  at  high  latitudes  and  two  stations  at  equatorial  latitudes).  The 
objective  of  this  study  is  to  investigate  and  implement  improvements  to  the  WBMOD  model 
based  on  more  extensive  data  sets  in  both  the  equatorial  and  high-latitude  regimes. 

During  the  first  two  years  of  the  study,  the  equatorial  section  of  the  \og{CkL)  model  in 
WBMOD  was  upgraded  using  data  sets  collected  under  the  auspices  of  the  Directorate  of 
Geophysics  of  the  USAF  Phillips  Laboratory  (PL)  [Secan  et  ai,  1993;  Secan  and  Bussey, 
1993].  This  report  summarizes  the  upgrade  to  the  high-latitude  section  of  the  WBMOD 
log(C'jtL)  model. 

2.  A  Brief  History  of  WBMOD 

This  document  marks  the  passing  of  the  name  WBMOD  for  this  model.  The  name  was 
appropriate  when  it  was  first  coined  in  the  early  1980s  (denoting  the  WideBand  MODel)  as 
the  model  was  largely  based  on  analysis  of  data  from  the  Defense  Nuclear  Agency  (DNA) 
Wideband  satellite  experiment  [Fremouw  et  ai,  1978].  With  inclusion  of  data  from  the 
subsequent  DNA  HiLat  and  Polar  BEAR  satellites  [Fremouw  et  al.,  1985;  Fremouw,  1986] 
and  from  the  Phillips  Laboratory  equatorial  scintillation  network  [Secan  et  ai,  1993],  the 
name  WBMOD  no  longer  holds  any  significance.  We  are  retaining  the  name  for  one  final 
version,  which  will  be  denoted  WBMOD  Version  13,  as  it  is  our  plan  to  develop  a  completely 
new  user  interface  for  the  model  as  part  of  a  new  computer  code  which  will  be  named 
SCINTMOD. 

As  a  matter  of  historical  record,  the  full  lineage  of  SCINTMOD  begins  with  the  BPMOD 
code  developed  at  SRI  International  in  the  early  1970s  [de  la  Beaujardiere  and  McNiel,  1971; 
Fremouw  and  Rino,  1973].  This  model  was  based  on  a  propagation  theory  from  Briggs  and 
Parkin  [1963]  and  a  global  model  of  {ANe)  (RMS  ANg)  derived  from  analysis  of  scintillation 
data  published  in  the  open  literature  [Fremouw  and  Bates,  1971].  The  successor  to  BPMOD 
was  a  model  called  lONSCNT  ]Rino  et  ai,  1977;  Fremouw  and  Rino,  1978]  which  replaced  the 
propagation  model  with  a  new  phase-screen  propagation  theory  [Rino  and  Fremouw,  1977] 
and  retained  the  BPMOD  model  for  {ANe).  The  first  version  of  WBMOD  was  developed 
in  the  early  1980s  [Fremouw  and  Lansinger,  1981],  which  replaced  the  propagation  model  in 
lONSCNT  with  the  completed  phase-screen  theory  as  described  in  Rino  [1977]  and  replaced 
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the  (AiVe)  model  with  a  model  of  the  path-integrated  irregularity  strength.  The  WBMOD 
model  has  passed  through  many  upgrades  and  improvements,  culminating  in  the  model 
described  in  this  document  and  in  Secan  and  Bussey,  [1993]. 

3.  High-Latitude  Model  Development 

With  completion  of  the  equatorial  region  of  the  model,  attention  has  been  moved  to 
upgrading  the  high-latitude  region.  A  very  different  set  of  data  will  be  used  in  this  part  of 
the  model-upgrade  effort  from  that  used  in  the  equatorial  upgrade.  As  with  the  equatorial 
database,  one  of  the  components  of  the  high-latitude  database  is  a  set  of  data  taken  from  the 
WIDEBAND  satellite  experiment,  collected  at  Poker  Flat,  Alaska.  The  remaining  data  were 
collected  from  the  HiLat  and  Polar  BEAR  satellites  from  four  ground  stations:  Bellevue, 
Washington;  Tromso,  Norway;  Ft.  Churchill,  Canada;  and  Sondre  Stromfjord,  Greenland. 
Table  lists  information  about  each  of  these  five  data  sets.  The  column  labeled  IPP  in  this 
table  is  the  geomagnetic  (apex)  latitude  at  the  350-km  ionospheric  penetration  point,  and 
is  the  average  for  each  data  set. 


Table  1:  Data  available  for  use  in  the  high-latitude  model  upgrade. 


Ground  Station 

Satellite 

IPP 

Dates 

Points 

Bellevue 

HiLat/PB 

54.1° 

1984-1987 

86,240 

Poker  Flat 

WIDEBAND 

65.9° 

1976-1979 

61,340 

Tromso 

HiLat/PB 

67.1° 

1984-1987 

512,490 

Ft.  Churchill 

HiLat/PB 

69.7° 

1984-1988 

Sondre  Stromfjord 

HiLat/PB 

74.0° 

1984-1988 

696,020 

Data  types  — 

S4  :  Intensity  scintillation  index  (RMS) 

(T0  :  Phase  scintillation  index  (RMS) 

T  :  Phase  PDS  spectral  strength  (IHz) 

p  :  Phase  PDS  spectral  slope 

The  basic  parameter  from  this  database  used  in  the  modeling  described  in  this  document 
is  CkL  calculated  from  the  slope  (p)  and  spectral  strength  (T)  parameters  derived  from  the 
phase  power-density  spectrum  (PDS).  The  method  used  to  calculate  CkL  from  T  and  p  was 
described  in  Section  3.2  of  Secan  and  Bussey,  [1993]. 
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3.1  Modeling  Philosophy.  We  were  able  to  use  a  somewhat  simpler  approach  to  modeling 
at  high  latitudes  than  that  used  for  the  equatorial  modeling  [Secan  et  al,  1993;  Secan  and 
Bussey,  1993]  because  the  distribution  of  iog(C'fcL)  in  the  high-latitude  region  is  very  nearly 
gaussian.  This  permits  us  to  model  the  average  log(C'jtT)  (the  PDF  parameter  fi)  and  the 
half-width  of  the  distribution  of  log(C'A:i)  (cr)  independent  of  one  another.  In  addition,  we 
have  access  to  the  full  \og{CkL)  PDF  since  we  are  using  T  and  p  to  derive  CkL,  whereas 
in  the  equatorial  modeling  we  used  SI{dB)  and  S4  to  derive  CkL,  both  of  which  lead  to 
distortions  of  the  high  end  of  the  log(C'fcL)  distribution  due  to  Fresnel  filtering.  Thus  our 
approach  to  the  high-latitude  modeling  is  to  develop  a  model  for  the  average  log(C'fcL)  to 
provide  the  PDF  parameter,  and  then  develop  a  second  model  for  cr. 

3.2  Variation  with  Geophysical  Conditions.  The  first  stage  in  the  log(C'A:T)  modeling 
was  to  (1)  establish  the  over-all  framework  of  the  high-latitude  log(C'itT)  model,  (2)  define  the 
features  to  be  included  in  this  model,  and  (3)  determine  the  variation  of  these  features  with 
local  time  and  with  variations  in  geophysical  conditions.  The  initial  plan  was  to  address 
these  items  separately,  but  it  became  obvious  early  on  that  the  complexity  of  the  high- 
latitude  behavior  would  require  us  to  deal  with  all  of  these  aspects  in  concert  rather  than  as 
independent  problems.  In  general,  we  moved  through  the  three  phases  of  this  stage  outlined 
at  the  start  of  this  paragraph  in  a  linear  fashion,  but  there  was  some  “doubling  back”  as  we 
found  aspects  of  the  framework  that  were  inadequate  or  additional  features  that  needed  to 
be  added. 

The  first  decision  concerning  the  framework  of  the  model  was  that  of  .the  coordinate 
“space”  in  which  the  model  would  be  defined.  As  described  in  Secan  and  Bussey  [1993], 
the  spatial  coordinate  system  is  the  distance  in  geomagnetic  latitude  from  the  equatorward 
boundary  of  auroral  electron  precipitation  (A)  and  geomagnetic  local  time  (r).  The  location 
of  the  equatorward  boundary  of  auroral  precipitation  is  calculated  from  a  model  of  that 
location  based  on  analysis  of  data  from  the  DMSP  SSJ  sensors  [Gussenhoven  et  al.,  1983]. 
This  model  can  use  either  an  observed  boundary  location  or  the  three-hour  Kp  index  from 
which  the  location  of  the  boundary  at  all  geomagnetic  local  times  is  produced.  We  found  early 
on  that  this  latitude  coordinate  ordered  the  scintillation  observations  well,  and  subsequent 
model  development  has  supported  this  finding. 

The  next  decision,  or  set  of  decisions,  was  the  definition  of  features  (boundaries,  peaks, 
depressions,  etc.)  to  be  included  in  the  new  high-latitude  model.  Our  decisions  were  based 
initially  on  visual  inspection  of  the  data,  decisions  that  were  revisited  and  modified  several 
times  during  the  modeling  process.  Figures  1,  2,  and  3  illustrate  the  types  of  displays  of  this 
large  data  set  that  were  used  in  this  definition  stage  and  in  later  modeling  work.  Figure  1  is 
a  gray-shade  plot  of  bin-averaged  log(C'jtT)  as  a  function  of  A  (radial  axis)  and  r  (azimuthal 
axis)  from  all  HiLat/Polar  BEAR  stations  for  all  SSN  and  all  Kp  values.  The  location  of  the 
precipitation  boundary  on  this  plot  (A  =  0)  is  indicated  by  the  solid  white  circle.  Figures  2 
and  3  show  these  same  data,  with  bin-averaged  log(C'fcL)  plotted  as  a  function  of  A  for  six 
geomagnetic  local  time  (GMLT)  periods. 
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Composite  -  All  Kp  /  All  SSN 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 


Figure  1:  Composite  shade  plot  of  log(C'A;L)  for  all  SSN  and  all  Kp.  The  radial  axis  is  the 
latitude  distance  from  the  SSJ/4  precipitation  boundary  (solid  white  line)  and  the  azimuthal 
coordinate  is  geomagnetic  (apex)  local  time.  The  solid  white  circle  indicates  the  zero  point 
for  the  radial  axis  {i.e.,  the  SSJ/4  boundary). 
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[Note:  This  and  all  other  gray-shade  plots  shown  in  this  document  are  somewhat  poor 
substitutes  for  the  corresponding  color  plots  actually  used  in  the  modeling  effort.  Some 
effort  was  made  to  find  a  good  gray-shade  range  that  would  highlight  features  in  these  plots 
as  clearly  as  they  are  in  a  color  plot,  but  many  features  that  are  very  prominent  in  color 
plots  do  not  stand  out  as  well  in  these  gray-shade  plots.] 

In  reviewing  these  figures  (and  as  a  result  of  iterations  during  the  modeling  that  followed), 
the  following  features  were  identified  and  included  in  the  new  high-latitude  log(C'fcL)  model: 


•  A  high-latitude  boundary  across  which  CkL  changes  by  over  an  order  of  magnitude. 
This  boundary  is  poleward  of  the  precipitation  boundary  by  a  distance  that  is  markedly 
larger  on  the  day  side  than  on  the  night  side. 

•  Enhancements  in  CkL  along  the  top  of  this  boundary  near  the  midnight  and  noon 
sectors.  (This  latitude  band  in  which  the  enhancements  occur  is  somewhat  loosely 
denoted  the  “auroral”  region  in  later  discussions.  While  it  is  clearly  within  the  au¬ 
roral  precipitation  region  in  the  midnight  sector,  it  is  significantly  poleward  of  such 
precipitation  in  the  noon  sector.) 

•  High  levels  of  CkL  in  the  noon-to-dusk  sector  of  the  polar-cap  region. 

•  A  depression  in  CkL  in  the  midnight-to-dawn  sector  of  the  polar-cap  region.  (This 
feature  is  denoted  the  polar  cap  “hole”  in  later  discussions.) 

•  A  transition  region  equatorward  of  the  sharp  transition  into  the  auroral  region  where 
the  CkL  levels  are  elevated  above  typical  mid-latitude  values  for  CkL.  This  region  is 
very  evident  in  the  0800  to  1200  GMLT  sector  in  Figure  2  and  the  1200  to  1600  GMLT 
sector  in  Figure  3.  (This  feature  is  denoted  the  sub-auroral,  or  SA  region,  in  later 
discussions.) 


After  several  iterations  through  the  process  of  determining  the  values  for  model  co¬ 
efficients  (described  later),  a  set  of  model  equations  was  constructed  that  reproduce  the 
variations  and  patterns  found  in  the  data.  The  basic  equation  for  calculating  log(C'fcL)  in 
the  high-latitude  region  is 


log(C,L)  =  C'o  +  (C'^-Co)exp 


/A-A^y 
\  W  ) 


(1) 


A  =  Xr 


Aj4 


(2) 


where  Ca  is  the  value  of  log(CA;T)  at  the  center  of  the  auroral  enhancement  region;  Co  is  the 
value  of  log(C'fcL)  either  in  the  polar  cap  region  (if  A  is  poleward  of  the  center  of  the  auroral 
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enhancement  region,  A^)  or  in  the  sub-auroral  (SA)  region  (if  A  is  equatorward  of  A^);  \m  is 
the  geomagnetic  latitude;  and  Aj4  is  the  geomagnetic  latitude  of  the  equatorward  boundary 
of  auroral  precipitation  as  specified  by  the  DMSP  SSJ/4  model  either  from  observations  of 
the  boundary  location  or  from  Kp.  The  center  latitude  of  the  auroral  enhancement  region  is 
given  by 


A.  = 


^A,N  +  ^A,D 


+ 


^A,N  —  ^A,D 


cos  (t-ta) 


(3) 


where  ta  is  defined  as  the  time  on  the  nightside  where  the  boundary  reaches  its  maximum 
equatorward  extent. 

The  Co  and  W  parameters  in  Eq.  (1)  are  calculated  as  follows: 

For  A  <  A^ 


Co  —  Cm  + 


2 


1  +  erf 


(4) 


and 


^A,SA,D  +  Wa,SA,N 


+ 


Wa,SA,D  —  ^A,SA,N 


COS  (t-ta) 


where  Cm  is  the  mid-latitude  log(C'fcL)  given  by 


(5) 


-  Cm  =  Cm,o  +  ^Cm,rR  +  ^Cm,t  cos(r ) 


(6) 


where  R  =  SSN /lOO  and  r  is  the  geomagnetic  local  time,  ACsa  is  the  difference  between 
the  mid-latitude  log(C'A;jC)  and  log(C'fcT)  in  the  sub-auroral  region,  given  by 


ACsa  =  ACsa,o  +  ACsa,t  cos  (r  -  tsa) 


(7) 


and  Asa  is  the  latitude  at  the  center  of  the  transition  from  the  mid-latitude  to  the  sub-auroral 
regions  given  by 


Asa  = 


Rsa,d  +  Asa,n 


+ 


RsA,D  —  Asa,N 


cos  (t- Tsa) 


(8) 
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where  tsa  is  defined  as  the  time  on  the  dayside  where  the  boundary  reaches  its  maximum 
equatorward  extent. 


Initially,  model  parameters  derived  in  an  earlier  effort  [Robins  et  al,  1986]  were  used  in 
the  mid-latitude  model  [Eq.  (6)].  In  deriving  the  parameters  for  the  high-latitude  model, 
however,  we  found  that  we  needed  to  adjust  these  parameters  in  order  to  obtain  the  best  fit 
to  the  observations.  This  is  discussed  further  in  the  next  section. 

For  A  > 


Co  =  Cpc 


Wa,PC,D  -\-Wa,PC,n'\  ,  \Wa,PC,D  —  Wa,PC,n'\  ,  N 

- +  - 2 - 


The  value  of  \og{CkL)  at  the  center  of  the  auroral  enhancement  region,  Ca,  is  calculated 
as  an  offset  from  the  log(C'fcZ')  value  in  the  polar  cap  region,  or 

Ca  =  Cpc  +  ACa  (11) 


where 


Cpc  =  Cpcfi  +  ACpc,H  X  G  (t ,Tpc,H,  Wpc  fj,  WpQ  jf^  x  F  (A,  Xpc,h) 


ACa  —  ACa,o  +  ACa,d  X  G  (^t,t^^p,Ata,d,^a,d^^a,d) 

+  ACa,N  X  G  (t,Ta,Nj  ATA,N,Wj^  j^,W^jsf^ 
The  two  functions  F  and  G  in  Eq.  (12)  and  Eq.  (13)  are  given  by 


FiX,Xo) 


if  X  <  a:o 
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Wf  =  (90  -  Xo)  ^(-ln|J|) 


(15) 


where  S  is  set  to  0.001,  and 


G{t,to,At,  W-,W+) 


'  1.0 

exp  — 

+  exp  [-  (^)"] 

T  within  of  To 

(16) 

otherwise 


where 


X+=(r-r„)-^ 
X-  ={t-  To)  +  ~ 


(17) 

(18) 


This  set  of  equations  then  defines  the  new  high-latitude  log(C'i;L)  model.  It  includes  29 
model  coefficients,  which  are  listed  in  Table  .  The  next  step  in  the  modeling  process  was  to 
establish  values  for  these  parameters,  and  to  determine  which  of  them  varied  as  functions  of 
either  SSN  or  Kp. 

The  method  used  to  determine  the  values  for  the  model  coefficients  is  similar  to  that  used 
in  the  Ionospheric  Scintillation  Analysis  Model  (program  IONS  AM)  developed  for  the  USAF 
Space  Forecast  Center  [Secan  et  al.,  1990].  The  model  is  treated  as  a  complex  function  that 
is  to  be  fitted  to  a  set  of  data  by  adjusting  all  or  some  subset  of  the  model  coefficients.  In 
the  lONSAM  program,  the  Levenberg-Marquardt  method  was  used  to  fit  the  high-latitude 
model  to  CkL  data  sets  derived  from  ground-based  scintillation  observations  and  in  situ 
density  measurements  from  the  DMSP  satellite.  Due  to  the  greater  complexity  of  the  new 
high-latitude  model,  it  was  decided  to  use  a  somewhat  simpler  method  for  fitting  the  model 
to  the  observations.  In  particular,  a  method  that  did  not  require  matrix  inversion  was  desired 
due  to  the  much  larger  number  of  model  coefficients  (a  number  that  increased  from  29  to 
51  when  variations  with  SSN  and  Kp  were  added  to  the  model).  The  method  selected  was 
a  simple  gradient  search  method  as  described  in  Chapter  11  of  Bevington  [1969]  (see  the 
discussion  of  program  GRADLS  on  pages  215-222).  This  method  was  implemented  in  such 
a  way  that  any  of  the  model  parameters  could  be  adjusted  in  the  fitting  process. 

In  order  to  speed  up  the  process  of  determining  the  model  coefficients,  the  fitting  pro¬ 
cedure  used  several  sets  of  bin-averaged  log(CfcL)  values  rather  than  “raw”  log(C'fcL)  values 
directly  from  the  four  station  databases.  A  total  of  thirteen  data  sets  were  generated:  one 
that  included  data  for  all  SSN  and  all  Kp  (shown  in  Figures  1  through  3),  six  for  all  Kp  for 
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Table  2:  Definition  of  the  coefficients  for  the  new  high-latitude  model. 


Symbol 


Description 


log(C'A;L)in  polar  cap  region  (base  term) 


^Cpc,H 

TpC,H 

^TPC,H 

^PC,H 

^PC,H 


^Ca,n 

^ta,n 

ta,n 

Wt^ 


Depth  of  polar-cap  “hole” 

MLT  of  the  center  of  the  polar-cap  “hole” 
MLT  width  of  polar-cap  “hole” 

MLT  half-width  of  polar-cap  “hole”  (-) 
MLT  half- width  of  polar-cap  “hole”  (-|-) 
Poleward  boundary  of  polar-cap  “hole” 


Alog(CfcZ/)  at  the  auroral  boundary 


Alog(CfcL)  at  dayside  maximum 
Width  of  dayside  maximum 
MLT  of  dayside  maximum 
MLT  half-width  of  dayside  maximum  (-) 
MLT  half-width  of  dayside  maximum  (-f) 


Alog(C'fcL)  at  nightside  maximum 
Width  of  nightside  maximum 
MLT  of  nightside  maximum 
MLT  half-width  of  nightside  maximum  (-) 
MLT  half-width  of  nightside  maximum  (-f-) 


^A,D  Dayside  extent  of  the  auroral  boundary 
Aa,n  Nightside  extent  of  the  auroral  boundary 
ta  MLT  of  maximum  nighttime  extent 

^A,SA,N  Half-width  of  SA  -4  A  transition  (nightside) 

Wa,sa,d  Half-width  of  SA  — >  A  transition  (dayside) 

^A,PC,N  Half-width  of  A  ^  PC  transition  (nightside) 

Wa,pc,d  Half-width  of  A  ->  PC  transition  (dayside) 


^CsA,o  Alog(C'A;L)  in  the  SA  region 
^Csa,t  Diurnal  variation  of  log(C'fcL)  in  the  SA  region 
Asa,n  Nightside  extent  of  the  SA  region 
Asa,d  Dayside  extent  of  the  SA  region 
tsa  MLT  of  maximum  dayside  extent 

WsA  Half-width  of  SA  transition 
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six  ranges  of  SSN,  and  six  for  all  SSN  for  six  ranges  of  Kp.  Grey-shade  plots  for  these  data 
sets  are  included  in  Appendix  A.  The  log(C'fcL)  values  in  these  data  sets  were  the  average 
log(C'A;L)  from  all  four  stations  in  A-r  bins  that  were  2°  wide  in  latitude  by  30  minutes  wide 
in  time.  Similar  binned  data  sets  of  SSN  and  Kp  were  generated  for  each  case  for  use  in 
determining  variations  with  SSN  and  Kp. 

Initial  values  for  the  30  model  coefficients  were  established  first  by  inspection  of  the  all- 
SSN /a\\-Kp  case,  then  by  fitting  the  model  to  this  case  with  all  model  parameters  included 
in  the  analysis.  Figures  4  and  5  show  the  results  of  this  initial  fit  in  format  similar  to  that 
used  in  Figures  2  and  3.  In  these  figures,  the  heavy  solid  curve  is  the  result  of  the  fitting 
procedure,  and  the  heavy  dashed  curve  shows  the  variation  of  log(C'fcT)  from  the  current 
WBMOD  high-latitude  model.  The  two  model  curves  are  not  as  smooth  as  one  would 
expect  from  the  smooth  functions  they  are  based  on  due  to  the  fact  that  the  bin-averaged 
SSN  and  Kp  values  valid  at  each  data  point  were  used  in  the  models  rather  than  values 
averaged  over  the  entire  set. 

With  the  base  values  established  for  the  model  coefficients,  the  next  step  was  to  ascertain 
which  of  the  coefficients  would  vary  with  geophysical  conditions  and  to  establish  relationships 
between  these  coefficients  and  the  appropriate  geophysical  indices.  We  first  addressed  the 
issue  of  variation  with  SSN.  The  fitting  procedure  was  executed  on  each  of  the  six  data  sets 
constructed  for  various  ranges  of  SSN  with  all  of  the  30  coefficients  included  in  the  fitting 
procedure.  The  variation  of  each  of  the  coefficients  as  a  function  of  SSN  was  examined, 
and  those  that  showed  no  clear  definable  variation  with  SSN  were  set  to  average  values  and 
excluded  from  future  adjustment.  The  fits  were  performed  again,  and  the  variation  of  the 
remaining  non-fixed  coefficients  were  examined.  After  several  iterations,  we  settled  on  4  of 
the  30  coefficients  that  clearly  (and  in  a  simple  fashion)  varied  with  SSN:  Cpcfl  and  ta,d 
varied  linearly  with  SSN ;  and  AC^.o  and  AGs^i.o  varied  as  SSN'^.  [Note:  The  model  actually 
uses  a  parameter  R  =  SSN / 100  for  the  SSN.] 

With  an  initial  model  for  the  variation  with  SSN  established,  the  same  procedure  was 
executed  using  the  six  data  sets  constructed  for  various  ranges  of  Kp.  In  this  analysis,  we 
found  that  12  of  the  30  coefficients  varied  with  Kpi  ACpc,H,  AC'a.o,  ACsa,  tpc,Hi  ta,n, 
Na,Di  ^a,n-,  Asa,d,  Wa,sa,n,  and  Wa,sa,d  varied  linearly  with  Kp\  and  Cpcfi,  ACa,n,  and 
ACa,d  varied  as  Kp^. 

[Note:  Care  was  taken  in  establishing  the  functional  forms  of  the  variation  of  model 
coefficients  with  SSN  and  Kp.  In  all  cases  where  the  variation  is  quadratic,  the  coefficient 
to  the  linear  term  is  opposite  that  of  the  coefficient  to  the  quadratic  term.  Thus,  in  all  of 
these  cases,  there  is  a  saturation  effect  in  that  there  is  a  threshold  SSN  (or  Kp),  defined  by 
the  minimum  or  maximum  of  the  curve,  above  or  below  which  the  coefficient  is  “frozen”  at 
the  value  obtained  at  the  threshold  SSN  (or  Kp).  This  precludes  the  model  from  “running 
wild”  for  abnormally  high  values  of  either  geophysical  parameter.] 

The  final  step  in  this  process  was  to  fit  the  model  simultaneously  to  both  the  SSN- 
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Composite  -  All  Kp  /  All  SSN 


Distance  from  Boundary 


Figure  4:  Variation  of  log(C'fcL)  with  the  latitude  distance  from  the  SSJ/4  precipitation 
boundary  for  six  geomagnetic  (apex)  local  time  ranges  for  all  SSN  and  all  Kp.  The  light 
solid  curve  indicates  the  observations,  the  heavy  solid  curve  the  new  high-latitude  model, 
and  the  heavy  dashed  curve  the  current  WBMOD  high-latitude  model. 
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invariant  data  sets  and  the  A^p-invariant  data  sets.  All  of  the  bin-averaged  log(C'fcA)  (and 
the  corresponding  bin-averaged  SSN  and  Kp)  from  the  12  data  sets  were  passed  to  the  fitting 
procedure  using  the  same  methods  described  in  the  last  two  paragraphs.  At  this  point,  the 
number  of  model  coefficients  had  increased  to  53  by  including  coefficients  to  the  linear  and 
quadratic  variation  of  some  of  the  model  coefficients  with  SSN  and  Kp.  After  a  number 
of  iterations,  it  was  found  that  the  variation  of  the  ACpc,H  coefficient  with  Kp  was  no 
longer  significant,  and  this  variation  was  dropped  from  the  model.  In  the  final  analysis,  the 
following  model  coefficients  were  found  to  be  functions  of  either  SSN,  Kp,  or  both: 


Cpc,o  —  f{R,Kp,Kp^) 
ACaa  =  f{R,R\Kp) 

aCa,d  =  /(/4,aV) 
ACa,n  =  f{Kp,Kp'^) 
ACsAfl  =  f{R.R\Kp) 

rpc,H  =  f(Kp) 
ta,d  =  f{R) 
rA,N  =  f{Kp) 

Ra,d  =  f{Np) 

Aa,n  =  f{Np) 

Rsa,d  =  f{Kp) 


Wa,sa,n  =  f{Np) 

Wa,sa,d  =  f{Kp) 

The  final  values  for  the  52  model  coefficients,  including  the  variations  with  SSN  and  Kp,  are 
listed  in  Table  .  Comparisons  of  the  resulting  model  and  the  current  WBMOD  model  with 
the  twelve  SSN  and  Kp  data  sets  in  the  same  format  as  Figures  4  and  5  are  presented  in 
Appendix  B. 

This  new  model  for  the  behavior  of  CkL  at  high  latitudes  is  quite  a  bit  more  complex 
than  the  old  model  (in  WBMOD  Version  12),  which  is  not  surprising  given  that  the  old 
model  was  based  on  analysis  of  data  from  a  single  site  (Poker  Flat,  Alaska)  in  only  two  local 
time  sectors  (pre-noon  and  pre-midnight).  The  variations  in  local  time  and  with  geophysical 
conditions  include  the  following: 


•  The  sub-auroral  (SA)  region,  which  may  be  associated  with  effects  from  the  plasma- 
pause  (although  without  a  more  in-depth  study,  this  connection  must  be  classified  as 
strictly  a  supposition),  is  almost  non-  existent  in  the  night  sector,  but  extends  up  to 
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Table  3:  Values  of  the  coefficients  for  the  new  high-latitude  model. 


Value 

Symbol 

Description 

31.32 

0.84 

Cpcfi 

log(C'fcT)in  polar  cap  region  (base  term) 
Variation  with  R 

-0.12 

Variation  with  Kp 

0.023 

Variation  with  Kp^ 

-0.50 

^Cpc,H 

Depth  of  polar-cap  “hole” 

1.80 

Tpc,H 

MLT  of  the  center  of  the  polar-cap  “hole” 

0.33 

Variation  with  Kp 

1.75 

^TPC,H 

MLT  width  of  polar-cap  “hole” 

3.46 

^PC,H 

MLT  half-width  of  polar-cap  “hole”  (— ) 

3.18 

76.69 

^PC,H 

MLT  half- width  of  polar-cap  “hole”  (-f) 
Poleward  boundary  of  polar-cap  “hole” 

IMiBM 

AC'a.o 

A\og{CkL)  at  the  auroral  boundary 

Variation  with  R 

0.12 

Variation  with  R? 

0.081 

Variation  with  Kp 

0.91 

0.037 

^Ca,d 

Alog(CfcL)  at  dayside  maximum 

Variation  with  Kp 

-0.027 

Variation  with  Kp^ 

0.31 

^ta,d 

MLT  width  of  dayside  maximum 

12.01 

ta,d 

MLT  of  dayside  maximum 

0.71 

Variation  with  R 

3.02 

MLT  half-width  of  dayside  maximum  (— ) 

3.69 

Wto 

MLT  half-width  of  dayside  maximum  (-f ) 

0.59 

0.20 

^Ca,n 

Alog(C'fcL)  at  nightside  maximum 

Variation  with  Kp 

-0.022 

Variation  with  Kp^ 

1.64 

^TA,N 

MLT  width  of  nightside  maximum 

1.31 

ta,n 

MLT  of  nightside  maximum 

0.20 

Variation  with  Kp 

2.91 

WX,N 

MLT  half-width  of  nightside  maximum  (— ) 

2.69 

wij. 

MLT  half-width  of  nightside  maximum  (-f) 
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Table  3:  Values  of  the  coefEcients  for  the  new  high-latitude  model  (continued). 


Value 

Symbol 

Description 

13.52 

^A,D 

Dayside  extent  of  the  auroral  boundary 

-0.025 

Variation  with  Kp 

2.10 

^A,N 

Nightside  extent  of  the  auroral  boundary 

0.94 

Variation  with  Kp 

23.37 

ta 

MLT  of  maximum  nighttime  extent 

5.65 

^A,SA,N 

Half-width  of  SA  A  transition  (nightside) 

0.051 

Variation  with  Kp 

6.01 

Wa,SA,D 

Half-width  of  SA  —¥  A  transition  (dayside) 

0.39 

Variation  with  Kp 

2.52 

^A,PC,N 

Half- width  of  A  PC  transition  (nightside) 

3.19 

^A,PC,D 

Half- width  of  A  — >•  PC  transition  (dayside) 

0.52 

^CsA 

Alog(C'jtT)  in  the  SA  region 

0.44 

Variation  with  R 

-0.11 

Variation  with 

0.20 

Variation  with  Kp 

0.19 

^Csa,t 

Diurnal  variation  of  log(C'ifeT)  in  the  SA  region 

-3.16 

^SA,N 

Nightside  extent  of  the  SA  region 

-16.72 

^SA,D 

Dayside  extent  of  the  SA  region 

1.74 

Variation  with  Kp 

13.69 

tsa 

MLT  of  maximum  dayside  extent 

3.86 

WsA 

Half-width  of  SA  transition 

15°  equatorward  of  the  precipitation  boundary  in  the  post-noon  sector.  The  scintilla¬ 
tion  level  in  this  region  is  a  function  of  both  SSN  and  Kp  (increasing  with  both),  and 
the  boundary  on  the  dayside  moves  toward  the  precipitation  boundary  with  increasing 

Kp. 

•  The  main  scintillation  boundary,  denoted  in  the  model  equations,  is  much  farther 
poleward  of  the  precipitation  boundary  on  the  dayside  than  it  is  on  the  nightside,  with 
the  peak  of  the  auroral  enhancement  occurring  roughly  2°  poleward  of  the  precipitation 
boundary  near  midnight  and  almost  14°  near  noon.  Both  the  location  of  the  auroral 
enhancement  and  the  half-width  of  this  feature  change  with  the  level  of  geomagnetic 
activity  {Kp),  broadening  on  both  the  dayside  and  the  nightside  and  moving  farther 
poleward  off  the  boundary  on  the  nightside  and  closer  to  the  boundary  on  the  dayside. 

•  The  auroral  enhancement  region  (that  latitude  band  at  the  “top”  of  the  main  increase 
in  CkL  with  latitude)  has  two  distinct  maxima:  one  in  the  midnight  sector  and  one 
in  the  noon  sector.  Both  maxima  increase  in  intensity  with  increasing  Kp,  while  the 
center  time  of  the  night  peak  moves  to  later  local  times  (starting  at  0121  GMLT)  with 
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increasing  and  the  time  of  the  day  peak  moves  to  later  times  (starting  at  1157 
GMLT)  with  increasing  SSN. 

•  The  polar  cap  region  shows  a  strong  enhancement  in  the  post-noon  sector,  and  a 
depression  (or  hole)  in  the  post-midnight  sector.  The  overall  CkL  level  in  the  polar 
cap  increases  with  increasing  SSN  and  decreases  with  increasing  Kp.  The  hole  region, 
which  may  be  associated  with  lower  levels  of  density  and  of  polar-cap  patches  due  to 
convection  dynamics  within  the  polar  cap  region,  is  located  at  0148  GMLT  and  moves 
to  later  local  times  with  increasing  Kp.  The  decrease  in  CkL  within  the  center  of 
this  region  does  not  appear  to  change  in  a  systematic  fashion  with  either  SSN  or  Kp., 
although  there  is  some  evidence  that  the  difference  between  the  average  CkL  in  the 
hole  region  and  the  average  CkL  in  the  rest  of  the  polar  cap  may  become  smaller  with 
increasing  Kp. 


3.3  Seasonal/Longitudinal  Variation.  The  present  high-latitude  model  for  log(Cjt-L)  in 
WBMOD  does  not  include  a  seasonal  variation,  as  no  seasonal  variation  was  found  in  the 
analysis  of  Wideband  data  collected  at  Poker  Flat  [Fremouw  and  Lansinger,  1981].  There 
was,  however,  strong  evidence  at  other  latitude  sectors  for  such  a  variation,  as  shown  by  Basu 
[1975]  and  Basu  and  Aarons  [1980]  (in  particular,  see  Figure  3  in  Basu  [1975]),  and  a  number 
of  studies  over  the  past  decade  have  focused  on  the  issue  of  seasonal  and  longitudinal  (or 
UT)  variations  in  the  polar  cap  [de  la  Beaujardiere  et  al.,  1985;  Buchau  et  al,  1985;  Sojka 
et  al.,  1994].  The  next  step  in  this  model  upgrade  was  to  look  for  evidence  of  a  seasonal 
variation  in  the  HiLat/Polar  BEAR  data  set  and,  if  found,  develop  a  model  for  the  variation 
to  be  included  in  the  upgrade. 

Such  a  variation  was  found  and  it  is,  for  the  most  part,  in  agreement  with  the  variations 
reported  by  Basu  [1975]  and  Basu  and  Aarons  [1980].  Figures  6  through  9  show  the  variation 
of  the  residual  \og{CkL),  defined  as  the  observed  log(CfeL)  minus  the  log(C'fcL)  derived  from 
the  model  described  in  the  previous  section,  as  a  function  of  the  time  elapsed  since  the  launch 
of  the  HiLat  satellite  for  each  of  the  high-latitude  HiLat/Polar  BEAR  stations  for  each  of  four 
geomagnetic  local  time  (GMLT)  sectors.  The  clearest  signals  in  these  figures  are  the  strong 
seasonal  variations  at  Ft.  Churchill  in  both  the  midnight  and  noon  sectors,  and  at  Sondre 
StromQord  in  the  noon  sector,  although  there  are  varying  levels  of  a  seasonal  variation  visible 
in  other  stations/time-sectors  as  well.  Note  also  that  the  variations  at  Ft.  Churchill  in  the 
noon  sector  are  opposite  in  phase  to  those  in  the  midnight  sector,  a  phenomenon  reported 
by  Basu  and  Aarons  [1980]  in  their  data  set  from  Narssarssuaq,  Greenland. 

As  an  initial  step  toward  including  these  variations  in  the  model,  each  of  the  twelve  curves 
shown  in  Figures  6  through  9  were  fit  using  a  multiple  regression  analysis  to  an  equation  of 
the  form 


Alog(CfeL) 


—  Co  A  cr 


(19) 
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log(CkL)  -  MODEL  log(CkL)  -  MODEL  log(CkL)  -  MODEL 


GMLT  Period:  2100  to  0300  (night) 


200  566  932  1298  1664  2030 

Mission  Day 


Figure  6:  Variation  of  the  residual  log((7feL)  as  a  function  of  days  since  the  launch  of  HiLat 
at  Tromso  (top  plot),  Sondre  Stromfjord  (middle  plot),  and  Ft.  Churchill  (bottom  plot)  for 
the  night  sector  (2100  to  0300  GMLT). 


log(CkL)  -  MODEL  log(CkL)  -  MODEL  log(CkL)  -  MODEL 


GMLT  Period:  0900  to  1500  (day) 


Figure  7:  Same  as  Figure  6  for  the  noon  sector  (0900  to  1500  GMLT). 
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log(CkL)  -  MODEL  log(CkL)  ^  MODEL  log(CkL)  -  MODEL 
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GMLT  Period:  0300  to  0900  (AM) 
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Figure  8:  Same  as  Figure  6  for  the  morning  (AM)  sector  (0300  to  0900  GMLT). 
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log(CkL)  -  MODEL  log(CkL)  -  MODEL  log(CkL)  -  MODEL 


GMLT  Period:  1500  to  2100  (PM) 


200  566  932  1298  1664  2030 


Mission  Day 


Figure  9:  Same  as  Figure  6  for  the  evening  (PM)  sector  (1500  to  2100  GMLT). 
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where  D  is  the  day  of  the  year.  The  SSN  dependence  was  included,  as  it  was  obvious  that  at 
least  a  few  of  the  stations  showed  a  modulation  of  the  seasonal  curve  as  a  function  of  SSN. 
(This  ^STV-dependence  will  be  discussed  further  at  the  end  of  this  section.)  The  results 
of  these  fits  are  shown  in  Figures  10  through  13  (the  fits  are  the  dotted  curves)  and  are 
summarized  in  Table  . 


Table  4:  Results  of  fitting  seasonal  variation  curves. 


Station 

Co 

cr 

CD 

Do 

Day  Sector  (0900  to  1500) 

TR 

ss 

FC 

-0.02 

-0.35 

-0.14 

0.10 

0.10 

0.57 

0.08 

0.22 

0.43 

76.1 

6.5 

14.8 

1.45 

0.68 

1.17 

PM  Sector  (1500  to  2100) 

TR 

SS 

FC 

-0.16 

0.05 

0.06 

0.16 

-0.02 

0.25 

0.14 

0.22 

0.05 

61.5 

20.8 

76.4 

1.15 

0.72 

1.13 

Night  Sector  (2100  to  0300) 

TR 

SS 

FC 

-0.17 

-0.12 

-0.15 

0.57 

0.31 

0.45 

0.13 

0.06 

0.29 

161.7 

167.7 
183.3 

0.57 

0.97 

0.45 

AM  Sector  (0300  to  0900) 

TR 

SS 

FC 

-0.09 

-0.07 

-0.12 

0.34 

0.00 

0.13 

0.14 

0.11 

0.07 

153.2 

138.6 

104.7 

1.15 

0.72 

1.13 

It  is  apparent  from  these  figures  that  there  is  definitely  a  seasonal  variation  in  high- 
latitude  scintillation  and  that  it  is  a  complex  function  of  local  time,  longitude  sector,  and 
latitude  regime.  We  have  chosen  to  focus  on  the  dominant  characteristics  rather  than  attempt 
to  model  all  of  the  subtleties  of  the  variations  seen  in  these  figures  for  two  reasons.  First,  two 
large  pieces  of  the  “puzzle”  are  missing:  we  have  no  observations  from  the  Russian  longitude 
sector,  and  only  limited  data  from  the  polar  cap  (Sondre  Stromflord  in  the  night  sector). 
Second,  we  have  found  it  difficult  to  reconcile  behavior  seen  at  various  stations  in  the  same 
time  and  longitude  sectors  within  the  context  of  the  present  modeling  effort. 

Figure  14  shows  the  distribution  of  stations  at  which  data  used  in  this  study  were  col- 
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GMLT  Period:  2100  to  0300  (night) 


200  566  932  1298  1664  2030 

Mission  Day 


Figure  10:  Data  from  Figure  6  (night  sector)  with  the  results  of  the  regression  analysis 
(dotted  curve). 
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GMLT  Period:  0900  to  1500  (day) 


200  566  932  1298  1664  2030 


Mission  Day 


Figure  11:  Data  from  Figure  7  (day  sector)  with  the  results  of  the  regression  analysis  (dotted 
curve). 
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log(CkL)  -  MODEL  log(CkL)  -  MODEL  log(CkL)  -  MODEL 


GMLT  Period:  0300  to  0900  (AM) 


Figure  12:  Data  from  Figure  8  (AM  sector)  with  the  results  of  the  regression  analysis  (dotted 
curve). 
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log(CkL)  -  MODEL  log(CkL)  -  MODEL  log(CkL)  -  MODEL 


GMLT  Period:  1500  to  2100  (PM) 


200  566  932  1298  1664  2030 


Mission  Day 

Figure  13:  Data  from  Figure  9  (PM  sector)  with  the  results  of  the  regression  analysis  (dotted 
curve). 
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lected.  The  longitudinal  “hole”  in  our  data  set  is  quite  obvious  here  —  we  have  no  observa¬ 
tions  covering  the  longitude  sector  from  45°E  to  180°E.  Our  main  question  is  whether  the 
variations  seen  at  Ft.  Churchill  will  also  be  seen  in  this  missing  sector,  or  if  there  will  be  no 
seasonal  variation  in  this  sector.  Based  on  works  cited  earlier  in  this  section,  we  have  chosen 
the  latter  option  of  limiting  the  seasonal  variation  to  the  US  sector.  This,  then,  leads  to  the 
next  problem:  at  what  longitudes  is  this  behavior  bound.  On  the  Altantic  side,  we  use  the 
fact  that  there  is  little  seasonal  variation  seen  in  the  Tromso  data  set,  and  have  set  the  east¬ 
ern  bound  to  half-way  between  Narssarssuaq  and  Tromso.  Initially,  we  were  inclined  to  use 
previous  analyses  that  showed  a  similar  lack  of  seasonal  behavior  at  Poker  Flat  to  use  that 
station  as  a  bound  on  the  Pacific  side.  First,  however,  we  generated  plots  similar  to  those 
shown  for  the  HiLat/Polar  BEAR  stations  for  the  Poker  Flat  (Wideband)  data  (night  and 
day  sectors),  which  are  shown  in  Figure  15.  The  variations  seen  in  this  figure  are  strikingly 
similar  to  those  seen  at  Ft.  Churchill.  This  presents  the  problem  of  deciding  where,  to  the 
west  of  Poker  Flat,  to  terminate  the  seasonal  variation  with  no  data  on  which  to  base  our 
decision.  For  now,  we  have  arbitrarily  set  the  western  bound  of  the  seasonal  variation  just 
to  the  west  of  Poker  Flat. 


Figure  14:  Location  of  stations  used  in  the  analysis. 

The  reconciliation  problems  involved  the  seasonal  variation  in  the  night  sector  as  seen  at 
Ft.  Churchill  and  Sondre  Strom:Qord  (Figure  10)  and  that  observed  at  Narssarssuaq  [Basu, 
1975,  Figure  3].  The  Narssarssuaq  observations  were  made  largely  in  the  auroral  region,  as 
were  the  bulk  of  the  observations  at  Ft.  Churchill.  The  Sondre  StromiQord  observations  (in 
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GMLT  Period:  2100  to  0300  (night) 


GMLT  Period:  0900  to  1500  (day) 


0  366  732  1098  1464  1830 


Wideband  Mission  Day 

Figure  15;  Variation  of  the  residual  log(C'feL)  as  a  function  of  days  since  the  launch  of 
Wideband  at  Poker  Flat  for  the  night  (top  plot)and  day  (bottom  plot)  sectors.  Also  shown 
is  the  regression  fit  to  these  data  (dotted  curves). 


the  night  sector)  were  made  largely  in  the  polar  cap,  as  were  a  subset  of  the  Ft.  Churchill 
observations.  Both  Ft.  Churchill  and  Narssarssuaq  show  strong  seasonal  behavior  in  the 
night  sector,  but  Sondre  Stromfjord  shows  little.  The  Ft.  Churchill  and  Sondre  Strom:Qord 
data  sets  were  further  subdivided  into  latitude  categories  to  see  if  the  seasonal  signature 
was  absent  in  the  Ft.  Churchill  observations  made  in  the  polar  cap  (designated  here  as 
observations  more  than  5°  poleward  of  the  auroral  enhancement  latitude  from  the  model) 
or  present  in  the  Sondre  Strom:Qord  observations  made  in  the  auroral  region  (less  than  5° 
poleward  of  the  auroral  enhancement  latitude).  The  polar-cap  subset  of  the  Ft.  Churchill 
data  continued  to  show  a  seasonal  variation,  while  a  very  weak  seasonal  variation  appeared 
in  the  auroral  subset  of  the  Sondre  Stromflord  data. 

After  trying  numerous  schemes  for  modeling  this  behavior,  we  have  settled  on  a  compro¬ 
mise  solution.  The  seasonal  variation  is  modeled  independently  for  the  polar  cap  region  and 
the  auroral  (and  sub-auroral)  region.  The  point  (in  latitude)  at  which  these  two  regions  are 
divided  is  defined  as  5.5°  poleward  of  the  auroral  enhancement  latitude,  with  a  transition 
width  set  to  damp  out  the  seasonal  variation  seen  at  Ft.  Churchill’s  latitude  just  poleward 
of  the  equatorward  edge  of  Sondre  StromQord’s  field-of-view.  Note  that  this  problem  was 
not  encountered  in  the  day  sector,  where  both  Ft.  Churchill  and  Sondre  Stromfjord  are 
equatorward  of  the  “auroral  enhancement”  and  showed  the  same  seasonal  variation.  The 
seasonal  variation  is  implemented  in  the  following  equation: 


/s(A,  4>mi  Tm:  D,  SSN)  =  fs,A  +  ifs,PC  “  fs,A)  -fl(A) 


(20) 


where  A  is  as  defined  in  Eq.  (2),  <j)m  is  the  geomagnetic  longitude  ,  is  the  geomagnetic 
local  time,  D  is  the  day  of  the  year,  fs,A  and  fs,pc  are  the  seasonal  variations  in  the  auroral 
and  polar-cap  regions,  respectively,  and  Fi  is  an  interpolation  function  given  by 


FM)  =  2 


(21) 


The  variation  in  the  aurora  region  is  given  by 


Ss,A  T  (^2  F2(Tto)] 
where  G  is  given  by  Eq.  (16),  F2  is  given  by 


1  -f  erf 


(Pm  Pm,s'\ 

V  j 


(22) 


(23) 
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and  the  Ci  and  C2  parameters  are  the  seasonal  variation  in  two  adjacent  local-time  sectors 
given  by 

SSN 

Cn  —  Cn,0  “I"  ('n,R  j^qq  ^  ^n,D  (-^  ^n,o)  5 

where  n  is  1  -  4  corresponding  to  the  four  GMLT  sectors  (1:  midnight,  2:  AM,  3:  noon,  4: 
PM).  The  Tm,s  times  in  Eq.  (23)  are  the  GMLT  values  that  define  the  boundaries  between 
these  time  sectors  (3.0,  9.0,  15.0,  and  21.0). 

Given  the  limited  amount  of  “pure”  polar-cap  observations  other  than  the  night  sector 
from  Sondre  Strom^ord,  we  have  set  the  fs,pc  term  in  Eq.  (20)  to  zero.  We  were  tempted 
to  put  in  an  ad  hoc  variation  based  on  recent  modeling  work  by  Sojka  et  al.  [1994],  but  have 
decided  to  leave  this  issue  to  be  addressed  at  a  later  time  using  a  data  set  collected  within 
the  polar  cap.  (Note:  A  data  set  that  would  be  useful  in  this  regard  is  the  PL/GPIA  Thule 
scintillation  data  set.) 

The  22  model  parameters  in  these  equations  were  derived  from  the  considerations  dis¬ 
cussed  earlier  in  this  section  and  from  the  behavior  at  Ft.  Churchill  and  Poker  Flat.  Table 
summarizes  the  values  of  these  parameters.  The  seasonal  term  calculated  from  Eq.  (20)  is 
then  added  to  the  model  value  for  log(C'fcL)  calculated  from  Eq.  (1). 

Earlier,  we  mentioned  that  an  5S7V-dependence  was  included  in  the  seasonal  variation 
as  the  residual-log (CfcT)  plots  showed  a  clear  variation  with  SSN,  in  spite  of  our  assumption 
that  we  would  be  removing  SSN  effects  by  removing  the  model  for  SSN  variation  described 
in  Section  3.2.  The  fact  that  the  seasonal  variation  does  depend  on  SSN  is  not  in  itself 
surprising,  particularly  if  this  effect  is  due  to  movement  of  solar-produced  plasma  in  and 
around  the  high-latitude  region.  Our  concern  is  that  our  analysis  presupposed  that  the 
seasonal  variation  was  independent  of,  or  separable  from,  the  SSN  variation  so  that  we  could 
analyze  for  them  separately  (a  technique  used  in  previous  WBMOD  modeling  work  [Fremouw 
and  Robins,  1985]).  Given  time  and  resources,  we  would  reopen  the  geophysical  variations 
analysis  and  approach  the  modeling  of  these  variations  together,  much  as  the  SSN  and  Kp 
variations  were  dealt  with  together.  As  neither  time  nor  additional  resources  are  available, 
we  will  implement  these  variations  as  they  are  described  herein,  with  a  caveat  that  there 
may  be  a  double  accounting  for  some  5S7V-dependent  variations. 

3.4  Modification  to  Mid-Latitude  log(C'fcZy)  Model.  As  stated  in  Section  3.2,  we  found 
that  we  needed  to  adjust  the  parameters  to  the  mid-latitude  log(C'fcT)  model  [Eq.  (6)]  in . 
order  to  obtain  the  best  fit  to  the  observations.  As  a  result  of  earlier  modeling  work,  the 
base  log(CfcL)  value  [parameter  Cm,o  in  Eq.  (6)]  was  set  to  30.24,  the  variation  with  R 
{ACm,r)  was  set  to  0.40,  and  the  diurnal  variation  (AC'm.t)  was  set  to  0.30  [Robins  et  al., 
1986].  The  new  values  for  these  parameters,  derived  by  including  the  coefficients  to  Eq.  (6) 
in  the  fitting  process,  are  29.01,  0.40,  and  0.37. 
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Table  5:  Values  of  the  coefficients  for  the  seasonal  variation. 


Value 

Symbol 

Description 

5.5 

As 

Seasonal  latitude  boundary  (auroral/polar  cap) 

1.0 

Ws 

Width  of  latitude  transition 

345.0 

Center  longitude,  seasonal  variation 

180.0 

Width  of  longitude  sector,  seasonal  variation 

30.0 

Width  of  longitude  transitions 

0.5 

■saB 

Width  of  local-time  transitions 

Co 

DC  term  (midnight  sector) 

DC  term  (AM  sector) 

-0.15 

DC  term  (noon  sector) 

0.05 

DC  term  (PM  sector) 

■E9 

cr 

SSN  term  (midnight  sector) 

mm 

SSN  term  (AM  sector) 

■il 

SSN  term  (noon  sector) 

0.25 

SSN  term  (PM  sector) 

0.40 

CD 

Seasonal  magnitude  (midnight  sector) 

0.05 

Seasonal  magnitude  (AM  sector) 

0.45 

Seasonal  magnitude  (noon  sector) 

0.05 

Seasonal  magnitude  (PM  sector) 

183.5 

Do 

Seasonal  phase  (midnight  sector) 

104.5 

Seasonal  phase  (AM  sector) 

15.0 

Seasonal  phase  (noon  sector) 

76.5 

Seasonal  phase  (PM  sector) 

The  reduction  in  the  lead  term  (C'm.o)  is  not  surprising.  When  the  initial  value  was 
derived,  it  was  derived  solely  from  Stanford  Wideband  observations  that  provided  a  higher 
value  for  this  parameter  than  did  the  Bellevue  HiLat  observations  then  available.  We  chose 
to  go  with  the  higher  Stanford-Wideband  value  because  we  only  had  a  small  amount  of 
Bellevue-HiLat  data  and  most  of  it  is  from  high  mid-latitude,  and,  if  the  model  were  to  be 
in  error,  we  wanted  it  to  err  on  the  high  side.  We  had  already  reduced  the  value  of  C'm.o  to 
29.32  in  Version  12  as  a  result  of  a  parametric  study  conducted  for  a  system  design  project, 
and  this  change  will  move  the  model  down  into  complete  agreement  with  the  (now  expanded) 
Bellevue-HiLat /Polar  BEAR  measurements. 

3,5  Variation  of  the  log(C'fcL)  PDF  Width.  The  value  for  the  half-width  of  the  log(C'fcL) 
PDF,  (T,  used  at  high  latitudes  in  WBMOD  Version  12  was  based  on  an  analysis  of  a  subset 
of  the  HiLat/Polar  BEAR  database  conducted  several  years  ago  [Secan  and  Fremouw,  1988]. 
Values  for  the  half-width  were  extracted  by  fitting  normal  distributions  to  discrete  PDFs 
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generated  from  the  various  data  sets.  It  was  found  that  a  value  of  1.0  for  cr  was  fairly 
representative  of  all  times  and  locations  in  the  data  set.  In  the  present  analysis,  values  for 
a  have  been  calculated  from  the  standard  deviation  of  log(C'fcL)  from  the  mean  log(C'fcZ/), 
which  we  will  denote  asD-  The  relationship  between  a  derived  in  the  previous  analysis  and 
Usd  is  Usd  =  cr/v^.  Thus,  the  current  value  for  u  of  1.0  corresponds  to  a  usd  value  of 
l/\^. 

The  results  of  the  present  study  do  not  appreciably  differ  from  those  obtained  in  the 
early  one.  We  have  looked  at  various  aspects  of  the  variation  of  usd  (with  latitude,  local 
time,  SSN^  and  Kp)  and  while  we  have  found  variations  away  from  a  mean  value  for  usdi  we 
have  not  found  any  systematic  variations  significant  enough  to  include  in  the  model.  The 
results  are  fairly  well  summarized  by  Figures  16  and  17.  The  variation  of  usd  as  a  function 
of  A  for  various  GMLT  sectors  is  shown  for  all  four  Hilat/Polar  BEAR  stations,  and  the 
present  (l/\/2,  dotted  line)  and  new  (0.8,  solid  line)  constant  values  for  usd  are  shown  for 
comparison.  The  only  change  we  intend  to  make  in  this  parameter  from  WBMOD  Version 
12  is  to  change  the  constant  value  from  1.0  to  1.13. 

3.6  Comparison  with  Narssarssuaq  Observations.  Figures  18  (from  Basu  [1975])  and 
19  provide  an  independent  check  of  the  new  high-latitude  model.  The  data  in  Figure  18  were 
taken  on  a  propagation  path  between  Narssarsuaq  and  ATS-3  during  the  years  1968  through 
1972  during  the  time  interval  2200  to  0200  GMLT  for  Kp  <  3°.  The  data  plotted  in  Figure 
19  were  generated  from  the  new  model  using  SSN  =  80  and  Kp  =  1.66.  The  agreement 
between  the  two  is  fairly  good,  with  the  most  obvious  differences  being  an  underestimate 
of  the  percent  of  time  that  S4  exceeded  30%  during  the  summer  months  and  the  tendency 
for  the  observed  seasonal  variation  to  be  less  symmetric  than  the  sine  function  used  in  the 
model.  Given  that  the  only  information  we  used  from  the  Narssarssuaq  observations  was  the 
fact  that  a  seasonal  variation  was  seen  at  Narssarssuaq,  we  are  very  pleased  with  the  level 
of  agreement  shown  in  these  two  figures. 

4.  WBMOD  Version  13 

The  new  models  for  high-latitude  scintillation  described  in  this  report  have  been  imple¬ 
mented  in  Version  13  of  WBMOD.  The  primary  differences  between  Version  13  and  Version 
12  are  the  replacement  of  the  high-latitude  model  for  log(C7A;Z')  with  the  model  described 
in  this  report  and  the  addition  of  the  DMSP  SSJ/4  precipitation  boundary  as  an  explicit 
optional  input.  The  User  Guide  written  for  Version  12  has  been  updated  to  reflect  the 
changes  made  in  version  13.  Copies  of  this  document  were  provided  to  PL/GPIA  and  to  the 
Air  Force  Space  Forecast  Center  (AFSFC).  The  AFSFC  also  received  a  copy  of  the  research 
version  of  WBMOD  Version  13  compiled  to  run  on  an  Intel  80386-class  computer. 

4.  Conclusion 

The  upgrade  to  the  equatorial  and  high-latitude  regions  of  the  WBMOD  CkL  model  has 
been  completed,  and  the  upgrade  has  been  implemented  in  the  research  version  of  WBMOD 
Version  13.  The  WBMOD  User’s  Guide  has  been  modified  to  include  changes  made  to  the 
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All  Stations  -  All  SSN  /  All  Kp 


-15.0  -10.0  -5.0  0.0  5.0  10.0  15.0  20.0  25.0 

Distance  from  Boundary 


Figure  16:  Variation  of  usd  with  the  latitude  distance  from  the  SSJ/4  precipitation  bound¬ 
ary  for  six  geomagnetic  (apex)  local  time  ranges  for  all  SSN  and  all  for  all  four  Hi- 
Lat/Polar  BEAR  stations  (solid  curve:  Bellevue;  dotted  curve:  Tromso;  short-dash  curve: 
Ft.  Churchill;  long-dash  curve:  Sondre  Stromflord).  The  heavy  dotted  line  indicates  the 
value  for  asD  corresponding  to  WBMOD  Version  12,  and  the  heavy  dotted  line  indicates  the 
new  value. 
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Figure  17:  Continuation  of  comparison  in  Figure  16. 
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Figure  18:  Observed  seasonal  variation  of  scintillation  at  Narssarssuaq,  Greenland  (from 
Basu  [1975]). 


Seasonal  Variation  at  Narssarssuaq 


Season 


Figure  19:  Modeled  seasonal  variation  at  Narssarssuaq.  The  solid  line  is  the  percent  of  time 
that  S4  >  30%  and  the  dotted  line  is  the  average  S^. 
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model  as  part  of  the  high-latitude  upgrade. 

While  this  concludes  the  planned  upgrades  to  the  WBMOD  model,  there  are  still  several 
open  questions  that,  if  resolved,  could  improve  the  model  further: 


•  How  does  log(C'fcL)  vary  with  Kp  and  SSN  in  the  polar  cap,  particularly  from  just 
poleward  of  the  cusp  region  to  poleward  of  the  field-of-view  at  Sondre  Stromfjord  in 
the  midnight  sector? 

•  What  is  the  seasonal  variation  in  the  Russian  sector? 

•  What  is  the  seasonal  variation  in  the  polar  cap? 


As  noted  in  the  section  describing  the  seasonal  variations,  there  is  a  data  set  that  could  be 
used  to  investigate  the  polar-cap  variation,  a  set  of  scintillation  observations  made  at  Thule, 
Greenland,  by  PL/GPIA  over  a  several-year  period.  (It  was  not  included  in  this  analysis 
due  to  limitation  of  resources.)  Modeling  of  the  seasonal  variation  in  the  Russian  sector  may 
need  to  be  based  on  analysis  of  satellite  data,  such  as  that  shown  in  Sojka  et  al.  [1994]  due 
to  a  lack  of  a  large  database  of  ground-based  scintillation  observations  from  that  sector. 
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Appendix  A.  Grey-Shade  Plots  of  High-Latitude  \og{CkL). 


The  figures  in  this  appendix  show  the  variation  of  \og{CkL)  from  the  composite  Hi- 
Lat/Polar  BEAR  database  as  a  function  of  the  latitude  distance  from  the  auroral  precipita¬ 
tion  boundary  (A)  and  geomagnetic  (apex)  local  time  (r).  The  location  of  the  precipitation 
boundary  (defined  as  A  =  0)  is  indicated  on  all  figures  by  a  solid  white  circle.  The  first 
figure  is  for  the  entire  data  set,  i.e.,  for  all  SSN  and  all  Kp.  The  following  six  figures  are  for 
all  Kp  for  six  ranges  of  SSN:  0  to  25,  25  to  50,  50  to  75,  75  to  100,  100  to  125,  and  >  125. 
The  final  six  figures  are  for  all  SSN  for  six  ranges  of  Kp:  Oo  to  lo,  lo  to  2o,  2o  to  3o,  3o  to 
4o,  4o  to  5o,  and  >  5o. 

These  figures  show  the  data  used  in  the  fitting  process  described  in  Section  3.2. 
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Composite  -  All  Kp  /  All  SSN 
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Figure  20:  Composite  shade  plot  of  log(C'fcL)  for  all  SSN  and  all  Kp.  The  radial  axis  is 
the  latitude  distance  from  the  auroral  precipitation  boundary  (solid  white  line)  and  the 
azimuthal  coordinate  is  geomagnetic  (apex)  local  time.  The  solid  white  circle  indicates  the 
zero  point  for  the  radial  axis  (i.e.,  the  auroral  precipitation  boundary). 
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Composite  —  SSN  Range:  0  to  25 
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Step  size:  2.00 


Figure  21:  Composite  shade  plot  of  \og{CkL)  for  all  Kp  and  0  <  SSN  <  25.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  —  SSN  Range:  25  to  50 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis;  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1 .00 

Step  size:  2.00 


Figure  22:  Composite  shade  plot  of  log(C'A;Z')  for  all  Kp  and  25  <  SSN  <  50.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  —  SSN  Range:  50  to  75 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 


Figure  23:  Composite  shade  plot  of  \og[CkL)  for  all  Kp  and  50  <  SSN  <  75.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  -  SSN  Range:  75  to  100 
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Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 


Figure  24:  Composite  shade  plot  of  log((7fcii)  for  all  Kp  and  75  <  SSN  <  100.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  --  SSN  Range:  100  to  125 
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Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 


Figure  25:  Composite  shade  plot  of  log(C'fcL)  for  all  Kp  and  100  <  SSN  <  125.  Axis 
definition  and  annotation  are  as  in  Figure  20. 
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Composite  --  SSN  Range:  >  125 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1 .00 

Step  size:  2.00 

Figure  26:  Composite  shade  plot  of  log(C'fc/i)  for  all  Kp  and  SSN  >  125.  Axis  definition  and 
annotation  are  as  in  Figure  20. 
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Composite  --  Kp  Range:  Oo  to  lo 
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Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 

Figure  27:  Composite  shade  plot  of  log(C'fcL)  for  all  SSN  and  Oo  <  Kp  <  lo.  Axis  definition 
and  annotation  are  as  in  Figure  20. 


49 


Composite  -  Kp  Range:  lo  to  2o 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1 .00 

Step  size:  2.00 

Figure  28:  Composite  shade  plot  of  log(Cfeii)  for  all  SSN  and  1°  <  Kp  <  2o.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  --  Kp  Range:  2o  to  3o 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 

Figure  29:  Composite  shade  plot  of  log(C'fcL)  for  all  SSN  and  2°  <  Kp  <  3o.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  --  Kp  Range:  3o  to  4o 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1 .00 

Step  size:  2.00 


Figure  30:  Composite  shade  plot  of  \og[CkL)  for  all  SSN  and  3o  <  Kp  <  4o.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  --  Kp  Range:  4o  to  5o 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 

Figure  31:  Composite  shade  plot  of  \og{CkL)  for  all  SSN  and  4o  <  Kp  <  5o.  Axis  definition 
and  annotation  are  as  in  Figure  20. 
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Composite  —  Kp  Range:  >  5o 


Latitude  axis:  55  +  (GMLAT-BNDLAT)  Time  axis:  Apex  Local  Time 

Range:  40.0  to  90.0  Step  size:  1.00 

Step  size:  2.00 

Figure  32:  Composite  shade  plot  of  log(C'fcL)  for  all  SSN  and  Kp  >  5°.  Axis  definition  and 
annotation  are  as  in  Figure  20. 
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Appendix  B.  Results  of  High-Latitude  log(C'jtL)  Modeling. 


The  figures  in  this  appendix  shows  the  results  of  the  fitting  process  described  in  Section 
3.2.  Each  pair  of  figures  corresponds  to  one  of  the  13  data  sets  shown  in  Appendix  A, 
and  are  presented  in  the  same  order  as  the  figures  in  Appendix  A.  These  figures  show  the 
variation  of  log(C'A;L)  as  a  function  of  the  distance  from  the  auroral  precipitation  boundary 
(A)  from  the  data  sets  (light  solid  curves),  from  the  new  high-latitude  model  (heavy  solid 
curves),  and  from  the  WBMOD  Version  12.02  high-latitude  model  (heavy  dashed  curves). 
This  variation  is  shown  for  six  geomagnetic  (apex)  local  time  sectors:  0000  to  0400  GMLT, 
0400  to  0800  GMLT,  0800  to  1200  GMLT,  1200  to  1600  GMLT,  1600  to  2000  GMLT,  and 
2000  to  2400  GMLT.  The  vertical  dotted  line  in  each  plot  indicates  the  location  of  the  auroral 
precipitation  boundary  (A  =  0). 
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Figure  33:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  \og{CkL)  (light  solid  line)  for  all 
SSN  and  all  /C. 
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Figure  34;  Continuation  of  comparison  in  Figure  33. 
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Figure  35:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  \og{CkL)  (light  solid  line)  for  all 
K-n  and  0  <  SSN  <  25. 
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Figure  37;  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcL)  (light  solid  line)  for  all 
Kr,  and  25  <  SSN  <  50. 
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Figure  38:  Continuation  of  comparison  in  Figure  37. 
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Figure  39:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  \og{CkL)  (light  solid  line)  for  all 
Kr,  and  50  <  SSN  <  75. 
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Figure  40:  Continuation  of  comparison  in  Figure  39. 
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Figure  41:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcL)  (light  solid  line)  for  all 
K„  and  75  <  SSN  <  100. 
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Figure  43:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcL)  (light  solid  line)  for  all 
K„  and  100  <  SSN  <  125. 
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Figure  44:  Continuation  of  comparison  in  Figure  43. 
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Figure  45:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  \og{CkL)  (light  solid  line)  for  all 
Kn  and  SSN  >  25. 
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Figure  47:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcL)  (light  solid  line)  for  all 
SSN  and  Oo  <  <  lo. 
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Figure  49:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcL)  (light  solid  line)  for  all 


SSN  and  lo  <  Kp  <  2°. 


72 


H 

o  32.0 
^  31,0 

U 


S  30.0 
0) 


29.0 

-15. 


5.0  10.0  15.0  20.0  25.0 

Distance  from  Boundary 

Figure  50:  Continuation  of  comparison  in  Figure  49. 
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Figure  51:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log{CkL)  (light  solid  line)  for  all 
SSN  and  2o  <  <  3o. 
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Figure  53:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcZ-)  (light  solid  line)  for  all 
SSN  and  3o  <  Kr,  <  4o. 
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Figure  54:  Continuation  of  comparison  in  Figure  53. 
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Figure  55:  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C'fcL)  (light  solid  line)  for  all 
SSN  and  4°  <  Kp  <  5°. 
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Figure  57;  Comparison  of  WBMOD  Version  12.02  (heavy  dashed  line)  and  SCINTMOD 
Version  1.0  (heavy  solid  line)  to  composite  observations  of  log(C’fcL)  (light  solid  line)  for  all 
SSN  and  K„  >  5o. 
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Figure  58:  Continuation  of  comparison  in  Figure  57. 
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